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ABSTRACT 



In this paper, a time integrated search for point sources of cosmic neutrinos is 
presented using the data collected from 2007 to 2010 by the ANTARES neutrino 
telescope. No statistically significant signal has been found and upper limits on 
the neutrino flux have been obtained. Assuming an E~ 2 spectrum, these flux 
limits are at 1 — 10xl0~ 8 GeV cm~ 2 s _1 for declinations ranging from —90° to 
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40°. Limits for specific models of RX J1713. 7-3946 and Vela X, which include 
information on the source morphology and spectrum, are also given. 

Subject headings: astroparticle physics; cosmic rays; neutrinos 



1. Introduction 



One of the main goals of the ANTARES telescope (lAgeron et al.ll20! ldl ) is the detection 
of cosmic neutrinos and the identification of their sources. Neutrinos only interact via the 
weak interaction and are stable, making them unique probes to st udy the high-energy uni- 
verse. The production of high-energy neutrinos has been proposed (IHalzen fc Hooperll2002bl ; 
Bednarek et al.ll2005bt ISteckerll2005bl ) for several kinds of astrophysical sources in which the 
acceleration of hadrons may occur. In the interaction of cosmic rays with matter or radia- 
tion, charged pions are produced. In the decay chain of pions, neutrinos are produced. The 
detection of neutrinos may give valuable information on the origin of cosmic rays. It would 
also settle the question of the hadronic versus lepto nic mechanism in seve ral sources from 
which high-energy gamma rays have been observed (IBerezhko et al.ll2008al ). 



The best neutrino flux upper limits up to PeV energies for the So uthern hemisphere have 
been e stablished by the ANTARES experiment using 2007-2008 data ([Adrian-Martinez et al. 
201 lbl ). In the present paper, this analysis is extended by adding two more years of data with 
the full configuration of twelve detection lines. Furthermore, the information on the amount 
of light produced in the events, which is a quantity correlated to the neutrino energy and 
which helps to distinguish the atmospheric neutrino background from a potential high-energy 
signal, is taken into account. 

The structure of this paper is as follows. In Section [2] the ANTARES detector is briefly 
described. Sections [3] and 0] describe the online selection and the simulation, respectively. 
The track reconstruction is explained in Section The selection of events is described 
in Section El Section [7] is devoted to the evaluation of the detector performance. The 
search method and the limit setting are described in Sections [8] and |9j Section [TO] shows 
how the search is improved by including the energy information. Results are presented in 
Section [TTJ A cross-check based on an alternative method is explained in Section [12j Finally, 
the conclusions are summarised in Section [T3"l 
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The ANTARES detector 



The operation principle of neutrino telescopes is based on the detection of the Cherenkov 
light induced by relativistic muons produced in the charged current (CC) weak interactions of 
high-energy neutrinos close or inside the detector. The information on the time and position 
of the detected photons is used to reconstruct the muon trajectory, which is correlated with 
the direction of the incoming neutrino. Other signatures are also possible, such as the 
cascades produced in the CC interactions of electron and tau neutrinos and in the neutral 
current interactions of all neutrino flavours. In this analysis muons induced by high-energy 
neutrinos are used. For these events, the detector acceptance is large due to the long muon 
range and the neutrino direction is derived with an accuracy of a fraction of a degree. 



The construction of the ANTARES detector (lAgeron et al.l l2011df) was completed in 
2008 , after se veral years of site exploration and detector RfcD (jAmram et al.l2002 ; Aguilar et al 



20051 . l2011eh . The detector is located at (42°48' N, 6° 10' E) at a depth of 2475 m, in 
the Mediterranean Sea, at 40 km from the French town of Toulon. It comprises a three- 
dimensional array of 885 optical mod ules (OMs) looking 4 5° downwards and distributed along 
12 vertical detection lines. An OM (jAmram et al.l 120021 ) consists of a ten-inch photomulti- 
plier (PMT) housed in a glass sphere together with its base, a special gel for optical coupling 
and a //-metal cage for magnetic shielding. The OMs are grouped in 25 triple ts (or storeys) 



on eac h line, except for one of the lines on which acoustic devices are installed ( lAguilar et al. 



2011a ) and which therefore contains only 20 optical storeys. The total length of each line is 
450 m, which are kept taut by a buoy located at the top of the line. The lower 100 m are 
not instrumented. The distance between triplets is 14.5 m and the separation between the 
lines ranges from 60 to 75 m. The lines are connected to a central junction box, which in 
turn is connected to shore via an electro-optical cable. Figure [T] shows a schematic view of 
the detector. 

The detector also includes several calibration systems. The lines slowly move due to the 
sea current (up to ~ 15 m at the top of the line in case of currents of 20 cm/s). A set of acous- 
tic devices together with tiltmeters and compasses along the lines are used to reconstruct the 



shape of the lines and orientation of the storeys every two minutes (lAdrian-Martmez et al. 



201ld ). The acoustic system provides the position of each optical module with a precision 
better than 15 cm. The time calibration is performed by means of a master clock on shore 
and a set of optical beacons (four along each line). This allows f or a calibration of th e time 
offsets of the photomultipliers with a precision better than 1 ns (lAguilar et al.ll201ld ). 



In this analysis, data from 29-1-2007 to 14-11-2010 are used. The total integrated live- 
time is 813 days, out of which 183 correspond to the five line period. Causes of loss of 
efficiency are the periods of high bioluminescence and sea operations. 
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Fig. 1. — Schematic view of the ANTARES detector, consisting of twelve mooring lines 
connected to the shore station through an electro-optical cable. 
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3. Online selection 

The charge and time information of all signals from the PMTs which exceed a pre-define d 



threshold voltage, typically the equivalent to 0.3 single photo-electrons (lAguilar et al.ll2007l ). 
are first digitised into 'hits' and then sent to shore where they are filtered by a farm of 
PCs. For this analysis, two different filter algorithms were used to select the events. Both 
are based on the assumption that the optical background processes such as potassium-40 
radioactive decays and bioluminescence are not correlated and induce single photo-electron 
hits. Hence, a first selection of the signal requires hits with a high charge (usually > 3 
photo-electrons) or coincident hits within a time window of 20 ns on separate OMs of the 
same storey (LI hits). The first trigger requires at least five LI hits compatible with a muon 
track in any direction. The second trigger is defined as the occurrence of at least two LI hits 
in three consecutive storeys within a specific time window. This time window is 100 ns in 
case that the two storeys are adjacent and 200 ns for next-to-adjacent storeys. In addition 
to the events selected by the trigger, the singles count rate of each OM is stored. 



Simulations 



Simulations are required for determining the acceptance and angular resolution of the 
detector, since in the absence of a source these quantities cannot readily be measured. 

The event simulati on starts with the generation of upgoi ng muon neutrino even ts using 
the GENHEN package <jBaileyH2002ah . which uses CTEQ6D JPumplin et al.l l2002ch parton 
density functions for computing the deep inelastic charged current neutrino scattering cross 
section. The events are weighted according to the cross section and their probability to 
survive the passage through the Earth. If the neutrino interaction occurs near the detector, 
the hadro nic shower res ulting from the break-up of the target nucleon is simulated using 
GEANT ( iBrunnerl 12003). Otherwise, only th e resulting muon is propagated to the detector 
using the MUSIC code (jAntonioli et al.lll997l ). Atmospheric muons reconstructed as upgoing 
are a source of background for a neutrino signal and their rejection is a crucial point in this 
analysis as will be described in Section El Downgoing atmospheric muons were simulated 
with the program MUPAGE (jCarminati et al.ll2008bl : iBecherini et al.ll2006bf ) which provides 
parametrised muons and muon bundles at the detector. 

Inside the detector, the Cherenkov photons emitted along a muon track and arriving on 
the OMs are simulated by sampling tabulated values of photon arrival times. The arrival 
time distributions have been derived by tracking individual photons t aking into account the 
measured absorption and scattering parameters (lAguilar et al.ll2005a[ ). 
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The PMT transit time spread (TTS) is simulated by a Gaussian smearing of photon 
arrival times. Optical backgrounds are added to the events according to the measured rates 
observed in the count rate data. Similarly, simulated hits from inactive OMs are deleted 
from the event. Sampling the count rate data from the runs selected for the analysis ensures 
that the simulation contains the same background and detector conditions as the analysed 
data set. The same trigger algorithms are applied to the simulation and the data. 



Aguilar et al. 


(2010). For the 


is considered ( 


3arr et al. 


2006a) 



5. Reconstruction 

Tracks a re recons t ructed from the hits in the triggered events using a multi-step al- 



gorithm (see iHeijboerl (120041 ) for a more detailed description). The initial steps provide a 
starting point for the final maximisation of the track likelihood. The likelihood is defined as 
the probability density of the observed hit time residuals, r, given the track parameters (po- 
sition at some arbitrary time and direction). The time residual r is defined as the difference 
between the observed and expected hit time for the assumed track parameters. 

It was found that the likelihood function has many local maxima and that the maximi- 
sation procedure needs to be started with track parameters close to the optimal values. The 
initial steps in the algorithm provide this near-optimal solution, estimating the track pa- 
rameters using increasingly refined score functions: a linear \ 2 fit, a so-called 'M-estimator' 
minimising g(r) = \/I + r 2 and a simplified version of the full likelihood fit. Each fit uses 
increasingly more inclusive hit selections based on the preceding stage. This sequence is 
started at nine different starting points to further increase the probability of finding the 
global optimum. 

The final likelihood function uses parametrisations for the probability density function 
(pdf) of the signal hit time residual, derived from simulations. The pdfs also include hits 
arriving late due to Cherenkov emission by secondary particles or light scattering. Further- 
more, the probability of a hit being due to background is accounted for. 

The quality of the track fit is quantified by the parameter 

A^-M^- + 0.1x(iV comp -l), (1) 
Nuts ~ 5 

which incorporates the maximum value of the likelihood, L, and the number of degrees of 
freedom of the fit, i.e. the number of hits, iVmtg, used in the fit minus the number of fit 
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parameters; N comp is the number of times the repeated initial steps of the reconstruction 
converged to the same result. In general, N comp = 1 for badly reconstructed events while it 
can be as large as 9 for well reconstructed events. 

The A varible can be used to reject badly reconstructed events, in particular atmospheric 
muons that are misreconstructed as upgoing. In addition, assuming that the likelihood 
function near the fitted maximum follows a multivariate Gaussian distribution, the error 
on the zenith and azimuth angles are estimated from the covariance matrix. The angular 
uncertainty on the muon track direction, (3, is obtained from these errors and can be used 
to further reject misreconstructed atmospheric muons as discussed in Section |6j 



6. Event selection 

Neutrino candidates are selected requiring tracks reconstructed as upgoing and applying 
selection criteria. These criteria were chosen following a 'blind' procedure, i.e. before per- 
forming the analysis on data. The effect of the selection cuts on data, expected background 
and signal efficiency are summarised in Table [TJ 

The estimated angular uncertainty on the muon track direction, (3, is required to be 
smaller than 1 degree. This cut rejects 47% of the atmospheric muons which are misrecon- 
structed as upgoing tracks. 

To further reject atmospheric muons that were misreconstructed as upgoing, the quality 
variable A is required to be larger than -5.2. This value is chosen to optimise the discovery 
potential, i.e. the neutrino flux needed to have a 50% chance of discovering the signal at the 
5<T significance level assuming an E~ 2 spectrum. 

Figure [2] shows the distribution of (3 for upgoing events with A > —5.2. The cumulative 
distribution of A for upgoing events is shown in Figure [3J The cut on the angular error 
estimate (3 is also applied. The excess of data compared to simulations at the lowest values 
of A is due to a non-simulated contribution of events consisting of solely optical background. 
Figure 0] shows the distribution of the reconstructed cosine of the zenith angle for both data 
and simulation. 

The final sample consists of 3058 neutrino candidate events out of a total of ~ 4 x 
10 8 triggered events. Simulations predict 358 ±179 atmospheric muons and 2408 ± 722 
atmospheric neutrinos, yielding a total expected events of 2766 ± 743. This is consistent 
with the observed rate within the quoted uncertainties (see Section H]). 
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Fig. 2. — Distribution of the estimate of the error on the direction of the reconstructed 
upgoing muon track after applying a cut on the quality variable A > —5.2. The red line shows 
the Monte Carlo atmospheric neutrinos, the purple line the Monte Carlo misreconstructed 
atmospheric muons and the black dots the data. The vertical dashed line with the arrow 
shows where the selection cut is applied (f3 < 1°). 
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Fig. 3. — Cumulative distribution of the reconstruction quality variable A for upgoing tracks 
which have an angular error estimate /3 < 1°. The bottom panel shows the ratio between 
data and simulations. The red line is for Monte Carlo atmospheric neutrinos, the purple 
line Monte Carlo misreconstructed atmospheric muons and the black dots the data. The 
vertical dashed line with the arrow shows where the selection cut is applied (A > —5.2). 
The purple and red bands show the systematic uncertainties on the simulations as explained 
in Section HI The green band in the bottom panel shows the total contribution of these 
uncertainties. 





data 


atm. \i atm. v E u 2 v 


Triggered events 


3.94 x 10 8 


3.06 x 10 8 1.54 x 10 4 100% 


Reco. upgoing events 


6.08 x 10 7 


2.98 x 10 7 1.24 x 10 4 61% 


Reco. upgoing events + (3 < 1° 


3.90 x 10 7 


1.57 x 10 7 8352 44% 


Reco. upgoing + j3 < 1° + A > -5.2 


3058 


358 2408 23% 



Table 1: Number of events before and after applying selection cuts for data (second column) 
and Monte Carlo simulations (third, fourth and fifth column). The last column shows the 
percentage of signal events assuming a neutrino flux proportional to an E~ 2 spectrum. 
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Fig. 4. — Distribution of the cosine of the zenith angle showing events with A > —5.2 
and {3 < 1°. The bottom panel shows the ratio between data and the total contribution of 
simulations. The simulated distributions are shown for atmospheric muons and neutrinos. 
Systematic uncertainties on Monte Carlo atmospheric muons and neutrinos are shown by 
the purple and red bands respectively. The green band corresponds to the sum of these 
uncertainties. The vertical dashed line with the arrow shows where the cut on the zenith 
angle is applied in order to select only upgoing events. 



13 



7. Detector performance 

The response of the detector to a neutrino signal proportional to an E~ 2 spectrum was 
obtained using the simulation described in Section @] and applying the analysis cuts. 

7.1. Angular resolution 

Figure [5] (left) shows the cumulative distribution of the angle between the direction of 
the reconstructed muon and that of the true neutrino. The median of this distribution is 
0.46 ± 0.10 degrees. Of the selected events, 83% are reconstructed better than 1 degree. For 
the data sample in which the detector was operational with all the 12 lines, the estimated 
angular resolution is 0.43 ± 0.10 degrees. The median of this angular error for the full data 
set considered in the analysis is shown in Figure (right) as a function of the true neutrino 
energy. 

The systematic uncertainty on the angular resolution quoted above has been estimated 
by varying the hit time resolution, At, in the simulation. Many possible effects can contribute 
to this resolution, including the PMT transit time spread, miscalibrations of the timing 
system and possible spatial misalignments of the detector. The hit time resolution directly 
impacts both the angular resolution and the number of events passing the quality criteria. 
Simulations using different At values were compared with data in order to determine the 
best agreement in the lambda distribution, and was obtained for At= 2.5 ns. This can be 
compared to the TTS of the PMT, which is 1.3 ns (standard deviation). However, the PMT 
time response is not Gaussian and the degraded resolution was found to partly account for 
the tails. A At of 2.5 ns yields the quoted angular resolution of 0.46° and is the value used 
in the plots shown in this paper. 

For At=3.4 ns, the simulations show a deterioration in angular resolution of 30% and 
the number of selected neutrino events in data exceeds the simulated neutrino signal by 2a, 
where a refers to the uncertainty on the atmospheric neutrino flux model. Hence, this value 
of At is excluded by the data. Assuming a linear dependency, this argument translates to a 
(la) uncertainty on the angular resolution of ~ 15%. 



The absolute orientation of the detector is known with an accuracy of ~ 0.1° ( lAdrian-Martmez et al 

2011ch . 
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Fig. 5.— (Top): Cumulative distribution of the angle between the reconstucted muon 
direction and the true neutrino direction for upgoing events of the whole data set. A neutrino 
spectrum proportional to E~ 2 is assumed. (Bottom): The median angle as a function of the 
neutrino energy E v . In these plots the cuts A > —5.2 and f3 < 1° are applied. 
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7.2. Acceptance 

The neutrino effective area, A„, is defined as the ratio between the neutrino event 
rate, R V (E V ), and the cosmic neutrino flux, $(E V ). The flux is assumed to consist of equal 
amounts of and v^. The neutrino effective area depends on the neutrino cross section, the 
absorption of neutrinos through the Earth and the muon detection (and selection) efficiency. 
It can be considered as the equivalent area of a 100% efficient detector. Figure |6] shows the 
effective function of the neutrino energy and declination, 8 U . 

The analysis is primarily concerned with cosmic sources emitting neutrinos with an E~ 2 
power law of the form 

dNv - - - fJ^GeV-WV 1 , (2) 



dE v dtdS T VGeVy 

where the constant <fi is the flux normalisation. The acceptance, A{5 U ), for such a flux, is 
defined as the constant of proportionality between <ft and the expected number of events in 
the source direction and can be expressed in terms of the effective 

A(6 v ) = 4>- 1 j JdtdE u Af(E u ,6 u )^- t . (3) 

The acceptance for this analysis is shown on Figure |6] (bottom) . For a source at a declination 
of -90(0)°, A = 8.8(4.8) x 10 7 GeV" 1 cm 2 s which means that a total of 8.8(4.8) neutrino 
candidates would be selected from a point-like source emitting a flux of 10~ 7 x (E u /GeV)~ 2 
GeV- 1 cm" 2 s _1 . 

To constrain the systematic uncertainty on the acceptance, a comparison between the 
atmospheric neutrino data and a simulation was performed, in which the efficiency of each of 
the OMs was reduced by 15%, which yields a 12% reduction of the signal events for an E~ 2 
flux. The atmospheric neutrino yield would instead be reduced by 40% to be compared to 
the 30% error on its flux normalisation. Therefore, the 15% uncertainty on the acceptance 
can be therefore considered a conservative choice. 



8. Search method 

Two alternative searches for point-like neutrino sources have been performed. The full- 
sky search looks for an excess of signal events over the atmospheric muon and neutrino 
background anywhere in the visible sky, i.e. in the declination range [-90°, +48°]. In the 
candidate list search, the presence of an excess of events is tested at the locations of the 51 
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sin(S) 



Fig. 6.— (Top): the neutrino effective area for the selected events as a function 
of the neutrino energy E v for three different declination bands. (Bottom): acceptance of 
the detector which is proportional to the number of events that would be detected and 
selected from a point-like source at a given declination assuming a flux of 10~ 7 x (E u /GeV)~ 2 
GeV~ 1 cm _2 s~ 1 as a function of the sine of the declination. 
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pre-defined candidate sources g iven in Table [2j They include the 24 source candidates from 



Adrian-Martinez et all (l2011bl ) and 27 new sources selected considering their gamma ray 
flux and their visibility at the ANTARES site as the selection criteria. Among the Galactic 
sources only TeV gamma ray emitters are considered. No such requirement is imposed for 
extrag alactic sources as Te V gamma ray s may be abso rbed by the Extragalactic Background 



Light (lNikishovlll962l ; lGould et al.lll966l : Ijellevlll966al ). 



The algorithm for the cluster search uses an unbinned maximum likelihood (IBarlow 



19901 ) which is defined as 

log£ s+b = Vlog[/i s x T(ipi{a s ,5 s )) x A/" s (A^ ts ) 

(4) 

where the sum is over the events; J 7 is a parametrisation of the point spread function, i.e. 
the probability density function of reconstructing event i at an angular distance ipi from the 
true source location (a s ,5 s ); B is a parametrisation of the background rate obtained from 
the distribution of the observed declination of the 3058 selected events; /i s and /i bg are the 
mean number of signal events and the total number of expected background events; N^ its is 
the number of hits used in the reconstruction. A/" s (A"^ its ) and A/" bg (A"^ its ) are the probabilities 
of measuring A^ its hits for signal and background respectively. The distribution of A^ its for 
data and Monte Carlo events is shown in Figure CD Figure M shows the distribution of N hits 
for signal function of the true neutrino energy. 

In the candidate list search, the sum in Equation H] incorporates the events located in a 
cluster within 20 degrees of the source. Events further away have T ~ and thus contribute 
a constant factor to the likelihood. In the full-sky search potentially significant clusters are 
first identified by selecting at least 4 events within a cone of 3 degrees diameter. Using a 
larger diameter or a bigger /lower number of required events increases the computation time 
without a significant improvement in the sensitivity. 

In the candidate list search the likelihood is maximised by numerically fitting the mean 
number of signal events, fi s , with the source location fixed. In the full-sky search the like- 
lihood maximisation yields the source coordinates and fi s for each cluster. After likelihood 
maximisation a test statistic, Q, is computed: 

g = log£™ b x -log£ b , (5) 

where log is the maximum value of the likelihood provided by the fit and log £ b is the 
likelihood computed for the background only case (/i s = 0). A large value of Q indicates a 
better compatibility with the signal hypothesis. In case of a full-sky search only the cluster 
with the largest value of the test statistic is considered. 
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Fig. 7. — Distribution of the number of hits used in the reconstruction, for the selected data 
(black dots), and the total Monte Carlo background contribution, i.e. atmospheric muons 
and atmospheric neutrinos (solid green line). The dashed blue line corresponds to the cosmic 
neutrino signal assuming an E~ 2 spectrum. The distribution is normalised to the integral of 
the total number of events. All the cuts described in Section [6] are applied. 
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Fig. 8.— Distribution of the mean of number of hits, (iVhits), as a function of E v for 
Monte Carlo signal events assuming a neutrino spectrum proportional to E~ 2 . The error 
bars represent the standard deviation of the N hits distribution. The final cuts of the analysis 
described in Section [6] are included. 
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9. Pseudo-experiment generation and limit setting 

In order to evaluate the sensitivity of the analysis, pseudo-experiments are generated 
simulating background and signal. Background events are randomly generated by sampling 
the declination from the parametrisation B. The right ascension is sampled from a uniform 
distribution. The simulation of the signal is performed by adding events around the coordi- 
nates of the source, sampling in this case the angular distance ipi of the event i to the source 
location from a three-dimensional distribution of the reconstruction error as a function of 
the declination and the number of hits. At this stage the systematic uncertainties on the 
angular resolution and on the absolute orientation of the detector are included by varying 
the simulated parameters of the events of each experiment, such as the zenith and azimuth 
angles, with the required uncertainty. An example of the distribution of the test statistic Q 
obtained by performing the search on a large number of pseudo-experiments for the full-sky 
search is shown in Figure [9] for the background only hypothesis and for experiments where 
several signal events are added to the background at a declination 5 = —70°. Figure ED also 
shows the values of Q corresponding to p-values of 2.7 x 10 -3 and 5.7 x 10~ 7 , i.e. 3a and 
5a. To compute the latter value, the distribution of the test statistic for the background 
only hypothesis has been extrapolated using an exponential fit. 

The median sensitivity and the flux upper limits at 90% confidence level (C L) are 



computed following the Feldman & Cousins prescription ( jFeldman fc Cousins! Il998h . The 



systematic uncertainty on the acceptance is taken into account in the computation. 



10. Discovery potential 

Figure [10] shows the probability of making a discovery at the 3a and 5a significance 
level, as a function of the mean number of signal events. The same curves are shown for 
a search which does not use N hits in the likelihood. The inclusion of the Nuts pdfs in the 
likelihood function reduces the number of events (and therefore the signal flux) needed for 
a discovery by ~ 25%. 

The worsening of the 3a and 5a discovery probability for a neutrino flux model with an 
exponential cut-off parametrised as dN/dE = <fi x (E u /GeV)~ 2 exp(—E u /E c ), with E c the 
cut-off energy, was estimated. In this case, for a source at a declination of 5 = —70°, the 
mean number of signal events needed for a 5a discovery assuming a cut-off energy E c = 1 
TeV is a factor 2 higher compared to that without an exponential cut-off. 

Simulations show that for a source with Gaussian extension cr sourco = 1° at a declination 
of 5 = —70°, the flux needed to claim a 5a discovery is a factor 1.2 higher compared to a 
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Fig. 9. — Distribution of the test statistic Q for the full-sky search. The full yellow histogram 
is for the background only experiments. The red, blue and green lines are for 3, 6 and 9 
signal events generated from a source with an E~ 2 spectrum at a declination of —70°. The 
vertical dotted lines show the values of Q corresponding to the 3a and 5a significance level. 
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mean number of signal events 



Fig. 10. — Probability for a 3a (red lines) and 5a (blue lines) full sky search discovery as 
a function of the mean number of signal events from a source at 5 = —70° with a neutrino 
spectrum proportional to E~ 2 . The dotted blue and red lines are for the likelihood described, 
the solid lines refer to the case where iVhits is not used. The horizontal dotted black line 
corresponds to the probability to make a discovery in 50% of the pseudo-experiments. 
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point- like source. 



11. Results 

A map in equatorial coordinates of the pre-trial significances of every point in the sky 
that is visible below the horizon at the ANTARES site is shown in Figure [TT1 




0> 



Fig. 11. — Sky map in equatorial coordinates showing the p- values obtained for the point- 
like clusters evaluated in the full-sky scan; the penalty factor accounting for the number of 
trials is not considered in this calculation. 

In the full-sky search the most significant cluster is located at (a, 5) = (—46.5°, —65.0°), 
where 5(9) events are within 1(3) degrees of this position. For this cluster the fit assigns 5.1 
as signal events, and the value of the test statistic is Q = 13.1. The corresponding p- value 
is obtained by comparing the value of the observed test statistic Q with the simulated Q 
distribution for the background only hypothesis. The post-trial p-value is 2.6%, which is 
equivalent to 2.2a (using the two-sided convention). 

The results from the search up the 51 a priori selected candidate sources are presented 
in Table [2] and shown in Figure [12j The most signal-like candidate source in the list is 
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HESS J1023-575. The maximum likelihood fit yield /i s = 2.0 and the test statistic value 
is Q = 2.4. The post-trial p-value of this cluster is 41%, compatible with a background 
fluctuation. Since no statistically significant cluster of events has been found, upper limits 
(Feldman & Cousins at 90% CL) for an E~ 2 flux are obtained for the candidate list sources. 
These limits and the corresponding 90% CL sensitivity are plotted in Figure [T2l as a function 
of the source declination. Also indicated are the published limits from other experiments. 



v 10"' 

(0 



E 

o 

> 
O 



CM > 

US 



10" : 



E 10 ' 



10" 



10° — 



• ANTARES (813 days) 
□ IceCube 40 (375.5 days) 
o MACRO (2300 days) 
Amanda-ll (1387 days) 



ANTARES sensitivity (813 days) 

- - IceCube 40 sensitivity (375.5 days) 
a Super-K (2623 days) 

o J 




-80 -60 



-40 -20 



20 40 60 80 
declination (degrees) 



Fig. 12. — Limits set on the E~ 2 flux for the 51 sources in the candidate list (see Table [2]). 
Upper limits, previously reported by other ne utrino experiments, on sources from both 
Northern and Sou thern sky are also included (lAmbrosio et al.l l200ll ; iThrane et al.l l2009al ; 
Abbasi et al.l 120091 ) . The ANTARES sensiti yity of this analysis is shown as a solid line and 
the IceCube 40 sensitivity as a dashed line (lAbbasi et al.l 1201 ll ). 



11.1. Limits for specific models 



Measurements of TeV gamma rays from the H.E.S.S. (lAharonian et al.l l2007al . 120061 ) 
telescope may indicate a possible hadronic scenario for the shell-type supernova remnant 
RX J1713. 7-3946 and the pulsar wind nebul a Vela X. The first observation of RX J1713.7- 
3946 with the Fermi Large Area Telescope ( lAbdo et al.l 1201 lal ) shows that the gamma-ray 
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Table 2. Results from the search for high-energy neutrinos from sources in the candidate 

list. The equatorial coordinates (a s , S s ) in degrees, p-value (p) probability and the 90% 
C.L. upper limit on the E~ 2 flux intensity 9O%CL in units of 10 _8 GeV _1 cm _2 s _1 are given 
(sorted in order of increasing p-value) for the 51 selected sources. 



Source name 


a.[°] 


*«["] 


P 


JjHU /oOL 


Source name 


«.["] 


Ss[°] 


P 




HESS J1023-575 


155.83 


-57.76 


0.41 


6.6 


SS 433 


-72.04 


4.98 




4.6 


3C 279 


-165.95 


-5.79 


0.48 


10.1 


HESS J1614-518 


-116.42 


-51.82 




2.0 


GX 339-4 


-104.30 


-48.79 


0.72 


5.8 


RX J1713.7-3946 


-101.75 


-39.75 




2.7 


Cir X-l 


-129.83 


-57.17 


0.79 


5.8 


3C454.3 


-16.50 


16.15 




5.5 


MGRO J1908+06 


-73.01 


6.27 


0.82 


10.1 


W28 


-89.57 


-23.34 




3.4 


ESO 139-G12 


-95.59 


-59.94 


0.94 


5.4 


HESS J0632+057 


98.24 


5.81 




4.6 


HESS J1356-645 


-151.00 


-64.50 


0.98 


5.1 


PKS 2155-304 


-30.28 


-30.22 




2.7 


PKS 0548-322 


87.67 


-32.27 


0.99 


7.1 


HESS J1741-302 


-94.75 


-30.20 




2.7 


HESS J1837-069 


-80.59 


-6.95 


0.99 


8.0 


Centaurus A 


-158.64 


-43.02 




2.1 


PKS 0454-234 


74.27 


-23.43 




7.0 


RX J0852.0-4622 


133.00 


-46.37 




1.5 


IceCube hotspot 


75.45 


-18.15 




7.0 


1ES 1101-232 


165.91 


-23.49 




2.8 


PKS 1454-354 


-135.64 


-35.67 




5.0 


Vela X 


128.75 


-45.60 




1.5 


RGB J0152+017 


28.17 


1.79 




6.3 


W51C 


-69.25 


14.19 




3.6 


Geminga 


98.31 


17.01 




7.3 


PKS 0426-380 


67.17 


-37.93 




1.4 


PSR B1259-63 


-164.30 


-63.83 




3.0 


LS 5039 


-83.44 


-14.83 




2.7 


PKS 2005-489 


-57.63 


-48.82 




2.8 


W44 


-75.96 


1.38 




3.1 


HESS J1616-508 


-116.03 


-50.97 




2.7 


RCW 86 


-139.32 


-62.48 




1.1 


HESS J1503-582 


-133.54 


-58.74 




2.8 


Crab 


83.63 


22.01 




4.1 


HESS J1632-478 


-111.96 


-47.82 




2.6 


HESS J1507-622 


-133.28 


-62.34 




1.1 


H 2356-309 


-0.22 


-30.63 




3.9 


1ES 0347-121 


57.35 


-11.99 




1.9 


MSH 15-52 


-131.47 


-59.16 




2.6 


VER J0648+152 


102.20 


15.27 




2.8 


Galactic Center 


-93.58 


-29.01 




3.8 


PKS 0537-441 


84.71 


-44.08 




1.3 


HESS J1303-631 


-164.23 


-63.20 




2.4 


HESS J1912+101 


-71.79 


10.15 




2.5 


HESS J1834-087 


-81.31 


-8.76 




4.3 


PKS 0235+164 


39.66 


16.61 




2.8 


PKS 1502+106 


-133.90 


10.52 




5.2 


IC443 


94.21 


22.51 




2.8 












PKS 0727-11 


112.58 


11.70 




1.9 
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emission seems to be co mpatible with a leptonic or i gin. Ho wever, composite models are also 
possible as discussed in I Zirakashvili &: Aharonianl (j2010al ). 



In iKappes et al.l (l2007bl ) the neutrino flux and signal rates are estimated for these ob- 
jects using the energy spectrum measured by H.E.S.S and by approximating the source 
extension with a Gaussian distribution. The spectrum for these models is shown in Fig- 
ure [13] Assuming these models, and taking into account the measured source extension, 
90% CL upper limits on t he flux norm alisation were computed for both sources. The Model 
Rejection Factor (MRF) (jHilll l2003aJ ) . i.e. the ratio between the 90% CL upper limit and 
the expected number of signal events, is also calculated. Figure [TBI summarises these results. 
For RX J1713. 7-3946 the upper limit is a factor 8.8 higher than the theoretical prediction. 
For Vela X the upper limit is a factor 9.1 higher than the model. In both cases these are the 
most restrictive limits for the emission models considered. 




10° 10 7 
(GeV) 




10 2 10 3 10 4 10 5 10 6 10 7 

E v (GeV) 



Fig. 13. — Neutrino flux models (dashed lines) and 90% CL upper limit (solid lines) for 
RX J1713.7-3946 (left) and Vela X (right). Also shown are the E~ 2 point -source limit (soli d 
horizontal blue lines) presented in Table [2j The models are taken from (IKappes et al.ll2007bl ). 



12. Cross-check with the Expectation-Maximisation method 



The results discussed in the previous section have been cross-checked using the E xpectation- 
Maximisation (E-M) algorithm applied to the problem of the search for point sources ([Dempster et al 
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19771 ; lAguilar et al.ll2008[ ). The E-M method is an iterative approach to maximum likelihood 



estimations of finite mixture problems, which are described by different probability density 
functions. In the case of a search for point sources the mixture model can be expressed as 
the sum of two components: 



log£ 



s-b = 5Z log [(~) x ^Oi(asA)) 
/it 



7V s (Ni 



hits/ 1 



+ (1 



/it 



(6) 



where the signal pdf (J 7 ) is modelled as a two-dimensional Gaussian, D is a polynomial 
parametrisation of the probability distribution of the events in declination; as in Equation 
HI the sum runs over all the events in the data set, /i t , and the number of hits is used to 
better discriminate between background-like and signal-like events. 

In comparison with the previous search method, the E-M algorithm uses a different 
likelihood description of the events and follows an analytical optimisation procedure that 
consists of two steps. In the expectation step the log-likelihood is evaluated using the current 
set of parameters describing the source properties. Then, during the maximisation step, a 
new set of parameters is computed maximising the expected log-likelihood. These parameters 
are the number of signal events attributed to the source, the source coordinates (in the full- 
sky search) and the standard deviations of the Gaussian describing the signal pdf. In this 
sense the E-M method has the freedom to adapt itself to the extension of the source. The 
test statistic used to determine the significance of a potential point source is obtained as in 
Equation [5j 



12.1. Results 



The most signal-like cluster found in the full-sky search is the same as that found by 
the search method described in Section |HJ The number of signal events estimated by the 
algorithm is /i s = 5.3. The observed value of the test statistic, Q = 12.8, or a larger one 
occurs in p = 2.6% of the background only experiments. No significant excess of events 
was found in the location of any of the candidate list sources. The lowest p-value is 0.87 
(post-trial corrected) and corresponds to the position of 3C-279. The results obtained with 
the two search methods described above are consistent. 
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13. Conclusions 

The results of a search for cosmic neutrino point sources with the ANTARES telescope 
using data taken in 2007-2010 have been presented. A likelihood ratio method has been 
used to search for clusters of neutrinos in the sky map. In addition to the position of the 
reconstructed events, the information of the number of hits has been used as an estimate 
of the neutrino energy. This improves the discrimination between the cosmic signal and the 
background of atmospheric neutrinos. 

Two searches have been performed: within a list of candidate sources and in the whole 
sky. No statistically significant excess has been found in either cases. In the full-sky search, 
the most signal-like cluster is at (a, 5) = (—46.5°, —65.0°). It consists of 9 events inside 
a 3 degrees cone, to which the likelihood fit assigns 5.1 signal events. The corresponding 
p-value is 0.026 with a significance of 2.2<r (two-sided). The most significant excess in the 
candidate list search corresponds to HESS J1023-575, with a post-trial p-value of 0.41. 90% 
CL Upper limits on the neutrino flux normalisation are set at 1-I0xl0~ 8 GeV cm -2 s" 1 in 
the range from 4 to 700 TeV (80% of the signal), assuming an energy spectrum of E~ 2 , and 
are the most restrictive ones for a large part of the Southern sky. These limits are a factor ~ 
2.7 better than those obtained in the previous ANTARES analysis based on the 2007-2008 
data. Limits for specific models of RX J1713. 7-3946 and Vela X, which include information 
on the source morphology and spectrum, were also given, resulting in a factor ~ 9 above the 
predicted fluxes. 
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